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ان اک و پرایت رصری یکی ازع ی بای ات کد مر ی سارت سای قر زان ای 
مختلف شده است. از طرفی جهت مبارزه با این پدیده و اجرای پروژه‌های حفاظت خاک نیاز به آگاهی از 
اطلاعات مکانی منابع رسوبات تولیدی در یک منطقه است. روش انگشت‌نگاری رسوب بر پایه ردیاب‌های 
ژئوشیمیایی. آلی. نسبت‌های ایزوتوپی و نیز استفاده از مدل‌های ترکیبی مختلف باعث شناخت سهم منابع مختلف 
رسوب در یک ناحیه می‌شود. در تحقیق حاضر با استفاده از ترکیب بهینه ردیاب‌های آلی. ژئوشیمیایی و نسبت‌های 
ایزوتوپی اقدام به تفکیک منابع مختلف فرسایش و تولید رسوب و سپس تعیین سهم این منابع با استفاده از مدل- 
های کالینز (011105)), کالینز اصلاح‌شله (رعصتاامی ۷ موتا (۷0/02 لاندور (100060۲) و اسلاتری 
)٩12167((‏ شد. جهت تعبین بهترین مدل در این‌باره از شاخص‌های 06۳ 16۷1517 و تطبیق نتایج با اندازه‌گیری 
واقعی مقادیر رسوب برپایه مدل 91۷1 و واحدهای شدت فرسایش استفاده شد. عناصر کرین, مس تیتانیوم و 
سیلیکون به‌عنوان ردیاب‌های تفکیک کننده واحدهای کاربری اراضی شناخته شدند و مدل ترکیبی کالینز اصلاح 
شده (عطنلا) 01 با شاخص‌های 600 ۹۹//٩‏ و 18/15۳ ۲//۰۷به‌عنوان بهترین مدل اتتخاب شد. نسایج 
اندازه‌گیری صحرایی میزان رسوب تولیدی از کاربری‌های مختلف اراضی نشان داد مراتع با پوشش‌های گیاهی 
روط و راهم نس ۲۱۸۷ ور ارا کت بای ی ای عاکه قراز خار ان مها کی او 
دادند بیشترین سهم در فرسایش و رسوب حوزه مربوط به مراتع ضعیف و تخریب يافته منطقه با درصد مساحت 
۵ و سهم ۵۷/۰۴ درصد می‌باشد و اراضی زراعی کمترین نقش را تولید رسوب منطقه دارند. همچنین نتایج 
مدل موتا با غبریب همبستگی ۰/۹۲۴ کمترین اخثلاف را با مقاذیر آندازه‌گیری‌شده داشت. 

کلیدواژه‌ها: انگشت‌نگاری, ردياب‌هاه رسوب. مدل‌های ترکیپی» تنگ بستانک. 
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۱۸ جغرافیا و مخاطرات محیطی شماره پیست و دوم 
۱. مقدمه 

لازمه اجرای برنامه‌های حفاظت خاک و کنترل رسوب. کسب اطلاع از اهمیّت‌نسبی منابع رسوب و سهم 
آن‌ها در تولید رسوب و در نتیجه مناطق بحرانی حوزه آبخیز است (کالینز و همکاران!۲۴۹:۲۰۱۰) پالازون و 
همکاران) ۲۰۱۵؛ ۲۰۶۸). کاهش تولید رسوب نیازمند اجرای روش‌های مناسب کنترل رسوب و حفاظت 
خاک در مناطق بحرانی منبع رسوب در حوزه آبخیز است (پاتریک و همکاران» ۲۰۱۵؛ ۲۱۱۸) و لازمدةٌ 
شناحت روش مناسب و اجرای صحیح برنامه‌های حفاظت‌خحاک شناخت و شناسایی اهمیّت‌نسبی منابع 
رسوب و میزان سهم هر یک از منابع رسوب (چن و همکاران۱۲۴:۲۰۱۶) در تولید رسوب است. بنابراین 
توسعه و ارتقاء روش‌های مدل‌سازی منشأیابی رسوب در مدیریت جامع حوزه‌های آبخیز ضروری است 
(نجفی و صادقی,۱۶۶:۱۳۹۲). در روش‌های سئتی با مشکلات نمونهبرداری مکانی و زمانی و تنگناهای 
اجرایی مواجه بودند و بعضی از آن‌ها به زمان و هزینه زیادی نیاز داشتند (کالینز و والینگ* ۱۶۱:۲۰۰۴). به- 
دلیل وجود مشکلات یاد شده در مورد روش‌های سنتی؛ روش‌های جدید باعنوان ردیابی» منشاء‌یابی و 
انگشت‌نگاری به‌عنوان روشی جایگزین و مناسب مورد توجه محققان مختلف قرار گرفته است (کویتر و 
همکاران,۲۰۱۳؛ ۲۵) و باعث شناخت و درک بهتری از فرسایش شده‌است (فرنز و همکاران۵۱۴:۲۰۱۴. 
انگشت‌نگاری رسوب براساس ارتباط بین میزان ردیاب‌های موجود در رسوبات و منابع متناظر رسوب در 
سطح حوزه استوار است که امکان تشخیص منابع بالقوه تولید رسوب و میزان سهم هر منبع در رسوب 
حروجی از حوزه براساس خحصوصیات انگشت‌نگاری را فراهم می‌کند. در این روش ویژگی‌های فیزیکی؛ 
ژئوشیمیایی و آلی رسوب و منابع رسوب برای تعیین منابع اصلی رسوب و اهمیت نسبی آنها مورد بهره- 
برداری قرار می‌گیرند (والینگ؛ ۱۶۰:۲۰۰۵). از جمله تحقیقات داخلی می‌توان به تحقیق نصرتی و همکاران 
(۱۳۹۴) اشاره کرد. نامبردگان از یازده عنصر به‌عنوان ردیاب‌های اوّلیه استفاده کردند که چهار عنصر منگنزه 
کربن. منیزیم و کلسیم به‌عنوان ترکیب بهینه جهت تفکیک سهم کاربری‌های مختلف اراضی در مدل ترکیبی 
استفاده شدند. نامبردگان سهم اراضی کشاورزی. مرتع و جنگل در تولید رسوب معلّق و بارکف را به‌ترتیب 
۴ ۰/۴ و ۱۶/۲ درصد اعلام کردند. غلامی و همکاران (۱۳۹۴) با استفاده از مدل ترکیبی به بررسی سهم 
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نامبردگان از دو عنصر کلسیم و منيزیم به‌عنوان ردیاب استفاده کردند. نامبردگان سهم کاربری‌های کشاورزی» 
مرتع و اراضی‌بایر را به‌ترتیب ۵۰/۵۲ ۸۲۶/۳۴ و ۲۳/۱۴ اعلام کردند. نجفی و صادقی (۱۳۹۲) با استفاده 
از تهیه نقشه سیمای فرسایش, انگشت‌نگاری و اندازه‌گیری میدانی سهم منابع مختلف تولید رسوب را در 
حوزة آبخیز ایده‌لو محاسبه کردند. نامبردگان نتایج اندازه گیری میدانی منطبق بر نتایج روش انگشت‌نگاری 
رسوب دانستند. قابل ذکر است که کلیّه تحقیقات داخلی از مدل ترکیبی کالینز (1۹۹۷) و یا کالینز اصلاح 
شده "(۲۰۱۰) استفاده کردند و تاکنون گزارشی مبنی بر استفاده از ساير مدل‌های ترکیبی همچون موتا 
(۲۰۰۴), لاندور (۲۰۱۰) و اسلاتری (۲۰۰۰) گزارش نشده است. در زمینه تحقیقات خارجی می‌توان به 
تحقیقات چن و همکاران (۲۰۱۶) بروی واحدهای کاربری‌اراضی چراگاه» مراتع با پوشش گیاهی متراکم» 
جنگل و خندق‌ها با استفاده از روش منشاءیابی رسوبات اشاره کرد. نامبردگان بیان کردند که جنگل‌های 
جوان (پوشش ضعیف) و چراگاه‌ها باسهم‌های ۸۵۰/۵ و ۲۵/۶ بیشترین سهم را در تولید رسوب در 
محدوده مورد مطالعه داشته‌اند. ونگ و همکاران (۲۰۱۴) با استفاده از رسوبات به دام افتاده از پشت سد در 
فلات لس چین توانستند تغییرات سهم منابع مختلف رسوب را با استفاده از رصد تغییرات کاربری‌اراضی 
پایش کنند. مرور منابع مختلف نشان می‌دهد که برای تعیین اهمیّت و سهم‌نسبی منابع رسوب از لایه‌های 
اطلاعاتی مختلی مانند کاربری اراضی, فرسایش‌های سطحی و زیرسطحی (راسل و همکاران"۲۰۰۱) و 
واحدهای سنگ‌شناسی (کالینز و والینگ, ۲۰۰۲) استفاده شده‌است. همچنین جمع‌بندی سوابق پژوهش مستند 
موجود دلالت بر عدم‌توجه کافی (بجز مطالعه صادقی و نجفی (۱۳۹۲)) سوابق یادشده به ارزیابی مقایسه‌ای 
نتایم حاصل از هر یک از شیوه‌های به‌کار گرفته‌شده و مقادیر اندازه‌گیری‌های صحرایی در سطوح حوزه 
آبخیز بوده که ضرورت انجام تحقیقات در این زمینه را تأیید می‌نماید. در تحقیق حاضر جهت مشحص کردن 
سهم و اهمیّت‌نسبی کاربری‌های اراضی موجود در منطقه تنگ‌بستانک استان فارس از مدل‌های ترکیبی کالینزه 
کالینز اصلاح‌شده. لاندور و اسلاتری و نیز اندازه‌گیری‌های صحرایی به‌همراه دو شاخص ارزیابی "06۳ و 
"07 استفاده شد تا با تعیین ترکیب مناسبی از ردیاب‌های قادر به تفکیک منابع رسوب کارآیسی روش 


انگشت‌نگاری ترکیبی» در واحدهای کاربری اراضی مدنظر قرار گیرد. 
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۳۰ جغرافیا و مخاطرات محیطی شماره بیست و دوم 
۲ منطقه مورد مطالعه 

منطقه مورد مطالعه در این تحقیق با مساحت "1610 ۸۱/۷ تحت عنوان حوزه آبخیز بهشت گمشده (تنگ 
بسقاتک) دن حلوف ۸ کیلمعری شمال‌طرب. فهرستان راز و در موق مق افیا ۳۳ ۱۵۲۰۳ ۷۶ 
۳ ۲ شرقی و ۳۳ ۱۶ ۳۰ تا ۱۸ ۲۵ ۳۰ شمالی واقع شده است. میانگین بارش سالانه این حوزه ۶۰۹ 
اتضازین ک‌کان: بختباری» کواترنر رازک و پابده گوریی می‌باشند. شکل | موقعیّت و راه‌های دسترسی به 


این منطقه را نشان می‌دهد. 
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شکل ۱ موقعیّت محدوده مورد مطالعه و راه‌های دسترسی به آن 


لقتنم استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... "۳ 
۳ مواد و روش‌ها 

۱ نمونه‌برداری 

در این مطالعه نقشه کاربری‌اراضی تحت کاربری‌های اراضی زراعی» جنگل, باغات و مراتع به‌ترتیب با 
مساحت‌های ۱۸۳۲/۶۱ ۰۴۶۷۶/۷۵ ۲۸۹/۹۵ و ۱۳۷۰/۶۸ هکتار به‌عنوان واحدهای مولّد رسوب و مطابق با 
روش الگوریتم حداکثر تشابه (1۷11) با تصویر ماهواره لندست ۸ سنجنده آمآل) (تیرماه۵۴) تهیّه شد. 
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در این مطالعه تعداد ۱٩‏ عنصر شامل ۱۵ عنصر ژئوشیمیایی و ۲ عنصر آلی و ۲ نسبت ایزوتوپی استفاده 
شد. عناصر ژئوشیمیایی شامل باریم کادمیوم کروم. مس لیتیم منگنز, نئودیوم نیکل فسفر سیلیکون, 


نگودیمیوم ۱۴۴/ نودیمیوم ۱۴۳ و ۲ ماده آلی نیتروژن و کربن بودند. 


انتخاب ترکیب بهینه ردیاب‌ها 

برای تعیین ترکیب ردیاب‌های منشأٌیاب از یک روش آماری دو مرحله‌ای استفاده شد. در مرحله اوّل در 
صورت عدمتبعیّت از توزیع نرمال» از آزمون غیر پارامتری آزمون ۲ یا آزمون کروسکال _ والیس استفاده 
گردید (پلازان و همکاران,۲۰۷۳:۲۰۱۵). در مرحله دوّم تعیین تابع تحلیل تشخیص و ضرائب آن برای کاهش 
خصوصیّات انتخابی اولیّه ردياب‌ها به طوری که دارای حداقل همبستگی و حداکثر توان تفکیک منابع رسوب 
باشند به‌کار گرفته شده است. بعد از محاسبه توابع تشخیص سطح معنی‌داری آماری آن‌ها با معیارهای 
2حصقا ۱۷115 و مربع فاصله ماهالانوبیس صورت گرفت. برای نشان دادن توان جداسازی یا طبقه‌بندی 
تحلیل تشخیص. معنی‌داری توابع؛ اختلاف گروه‌ها و ماتریس طبقه‌بندی بررسی شده است. درصد طبقه‌بندی 
صحیح مهمترین معیار بررسی توان طبقه‌بندی تحلیل تشخیص می‌باشد که برای تعیین آن از ماتریس طبقه- 
بندی استفاده شده‌است. در ضمن به منظور ارزیابی نتایج تحلیل تشخیص از روش ارزیابی متقابل استفاده به 
عمل آمده است. ترکیب مناسب به دست آمده از ردیاب‌ها در این مرحله برای تعیین سهم منابع رسوب در 
مرحله بعد استفاده می‌شود (والینگ و همکاران»۲۰۰۸؛ ۲۷۹). یکی از پیش‌فرض‌های روش تحلیل تشخیص. 
نبود هم‌خطی چندگانه قوی بین متغیُرها می‌باشد. برای تشخیص هم‌خطی چندگانه از عدد تحمّل و عکس آن 
(عامل تورم واریانس) استفاده شده است (عباسی و همکاران ۱۳۸۹؛۳۷). عدد تحمّل برابر یا کوچکتر از ۰/۱ 
و عامل تورم واریانس برابر با بزرگتر از ۱۰ بیانگر هم‌حطی قوی بین متغیّرهاست و برای مقابله با مشکلات 
هم‌خطی چندگانه قوی» ترکیب‌هایی از عناصر طوری انتخاب شدند که کمترین هم‌خطی در بین آن‌ها (عامل 
تورم واریانس حداکثر ۱۰) وجود داشته باشد. 
۲ ۳ مدل‌های ترکیبی و کمّی 

تعین سهم منابع رسوب بر اساس مقایسه خصوصیّات نمونه‌های رسوب و منابع رسوب با استفاده از 
روش‌های کمّی با به‌عبارتی مدل‌های چندمتغیره ترکیبی صورت گرفت. برای حل این مدل‌ها برنامه‌نویسی 
خطی برای حل و بهینه کردن مقادیر سهم‌نسبی استفاده می‌شود (حددادچجی و همکاران,۲۰۱۳؛ ۸۵۶٩‏ والینگ و 
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۳ استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... ۳۳ 


شود: 
رابطه (۱) 


که در آن ز) برابر با مقدار برآوردی ردیاب ام برای نمونه‌های رسوب (10...» ۲ و 1-۱ زنة مقدار 
میانگین ردیاب لام در واحد کاربری اراضی لام (۲...0 و ۱-ز) ز میانگین مشارکت‌نسبی واحد کاربری- 
اراضی [ام (که نامعلوم بوده و به‌وسیله مدل برآورد می‌شود)؛ 9 تعداد واحدهای کاربری اراضی و ۲ تعداد 
ردیاب‌های انتخاب شده (ترکیب بهینه) می‌باشد. از آن جایی که فرض می‌شود ردیاب‌ها و ویژگی‌های آنان 
در نمونه‌های رسوب کاملاً وابسته به مواد منبع واحدها می‌باشند و این ردیاب‌ها تنها از این واحدها ناشی 


می‌شوند. بنابراین نتایج حاصل. مشروط به دو محدودیت (ز۳) زیر می‌باشند: 


رابطه (۲) 1 > ,۲ > 0 
رابطه (۳) و وحزا 


به این علّت که در چنین مدل‌هایی, رابطه (۱) برای تمام ردیاب‌ها تکرار می‌شود بنابراین تعداد معادلات 
بیشتری نسبت به متغیرهای مجهول وجود خواهد داشت. برای حل این معادلات از روش بهینه‌سازی تابع 
هدف و تکرار عملیّات استفاده می‌شود که در این پژوهش نیز از کمینه‌سازی مجموع مربعات باقی‌مانده با 
روش لوکال ([1,06) به‌عنوان تابع هدف استفاده شد. قبل از تعیین سهم باید اصلاحاتی از جمله ضریب 
وزنی ویژه هریک از ردیاب‌ها (1) که هریک از خصوصیّات منشأٌیاب را با چند تکرار در آزمایشگاه 
اندازه‌گیری کرده و با محاسبه عکس ریشه دوم واریانس داده‌های استاندارد شده به دست می‌آید» را محاسبه 
کرد و نیز ضریب کربن آلی (:60) (برای برطرف کردن مسئله غنی شدن يا تهی شدن رسوب) که با تقسیم 
میزان کربن آلی رسوب به متوسط میزان کربن آلی هریک از منابع رسوب بدست می‌آید رامحاسبه کرد و 
سپس اصلاحات را به مدل چندمتغیرةُ ترکیبی وارد کرد. همچنین برای تعیین ضریب تصحیح اندازه ذرات 
(2) با توجه به طبیعت انتخابی فرآیندهای فرسایش و انتقال رسوب باعث غنی شدن یا تهی شدن رسوب از 
نظر بعضی از اندازه‌ها می‌شود. در اثر این فرآیند. توزیع اندازه ذرات نمونه‌های رسوب و خحاک منشاً متفاوت 
خواهد بود؛ بنابراین مقایسه بین حصوصیّات خاک منابع رسوب و نمونه‌های رسوب که از نظر اندازه ذرات 
اختلاف دارند. صحیح نخواهد بود. جهت رفع این عیب و فراهم آوردن امکان مقایسه مستقیم بین نمونه‌های 
رسوب و نمونه‌های خاک باید از روش‌های اصلاحی استفاده شود. در این تحقیق بااین فرض که 
خصوصیّات منشأیاب به یک اندازه تحت تأثیر فرآیند غنی شدن يا تهی شدن قرار می‌گیرند و به دلیل نقش 
زیاد پارامتر سطح‌ویژه در میزان خصوصیات منشاأیاب (ردیاب‌ها) و انعکاس ترکیب اندازه ذرات» از نسبت 


۳۴ جغرافیا و مخاطرات محیطی شماره پیست و دوم 
سطح‌ویژه هر یک ازنمونه‌های رسوب به متوسط سطح ویژه نمونه‌های خاک هر یک از منابع رسوب برای 
اصلاح خحصوصیّات منشأٌیاب منابع رسوب استفاده شد. با حل مدل‌های زیر سهم تمام واحدهای کاربری- 
اراضی برای هر نمونه رسوب به‌صورت جداگانه به‌دست آمد که در نهایت با میانگین گیری از سهم واحدها؛ 


مقدار متوسط سهم هر واحد کاربری اراضی در تولید رسوب تعیین شد. 


رابطه (۴): (حدادچی و همکاران/۵۷۰:۲۰۱۳) 8 - وزنگ وح زره ]ور < ٩‏ 
رابطه (۵) ت(ربکر ی ۳۳ -.0), - ۳ 
4 

رابطه (ع) / رک و شورس سر گم _ 

3 / 

رعه تست زد | ۳ 

رابطه (۷): (پاتریک,۲۱۲۱+۲۰۱۵) ۹ 1 سس سح خرب2) سک -» ۳ 
رابطه (۸): (کالینز و همکاران۲۰۱۰؛ ۲۵۵) و ی عسک -, 


در روابط بالا ز:۷۸6 واریانس ویژگی ردیاب [ام منبع رسوب آام. 1۶ مجموع مربعات باقیمانده. ز) میزان 
اندازه‌گیری‌شده ویژگی (ردیاب) ام در نمونه رسوب. :۷۷ ضریب وزنی ویژه ردیاب لام. 2 ضریب اصلاحی 
اندازه ذرات برای فرسایش (منبع) [ام ز) ضریب تصحیح ماده آلی برای فرسایش (منبع) [ام» زنگ میزان 
ویژگی (ردیاب) لام در منبع [ام 0 تعداد منابع» ظ تعداد ردیاب‌ها. 130 تعداد نمونه‌ها در هر منبع» ز26 سهم 
منبع در فرسایش و تولید رسوب نز9۷٩‏ وزن درون منبع [برای مقادیر مختلف ردیاب‌های منشآیاب 1 برای 
تست دقت بهینه‌سازی مقدار *1 (مجموع مربعات باقیمانده) از روش پیشنهادی کالینز و همکاران (۲۰۱۲) 


مطابق با رابطه ٩‏ و برآورد خطای نسبی مدل مطابق با رابطه ۱۰ استفاده شد. 


رابطه ۹( 9( خ بت 
۱ رو -1 - 60 
رابطه (۱۰): (درپریش۵۲:۱۳۹۴) :2 - :0)ی 2 _ - 5 


نت( کِ )کب( 


در این روابط 1 تعداد خحصوصیّات ردیاب» 9 مقدار اندازه گیری‌شده ردیاب ام در نموبه رسوب» زز5 


مقدار فیان‌کین ردیاب لام در کاربری [ام. 2 سهم کاربری [ام. مصل) قیال کین مقادیر مشاهده‌شده 0 مقادیر 


1 ۳۱۵۸002021 6 ۰ 


بال‌شتیم استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... ۲۵ 


مشاهده‌شده و :0 مقدار برآوردی حاصل از مدل‌هاست. هرچه میزان شاخحص 20(۳) بیشتر باشد و مقدار 
۷7 کمتر باشد» دفت مدل بیشتر است. 


اندازه‌گیری‌های صحرایی 

به‌منظور ارزیابی نتایج به‌دست آمده از کاربرد مدل چندمتغیّره ترکیبی لازم بود تا مشاهدات صحرایی و 
اندازه‌گیری‌های مستقیم صحرایی نیز صورت پذيرد. بر همین اساس. با استفاده از روش دفتر مدیریت اراضی 
(1/۷1) اقدام به تهیه واحدهای کاری نیمه‌همگن شد. با بررسی‌های صحرایی و اصلاح و به‌روز رسانی 
کاربری‌ها و واحدهای نیمه‌همگن واحدهای همگن به‌دست آمدند. در راستای افزایش دقت نقشه ایجاد شده 
از روش ارائه شده توسط صادقی (۲۰۰۵) استفاده‌شد به‌طوری که طبق امتیازات حاصل از روش ]1,۷ طبقه 
فرسایشی برای هر واحد از مقاوم به فرسایش (امتیاز ۰ تا۲۰) فرسایش خیلی‌شدید (امتیاز ۸۰ تا ۱۰۰) طبقه- 
بندی شد و به‌ترتیب براساس شلات فرسایش موجود در واحد کاری حروف ۸ تال احتصاص داده‌شد که این 
حروف شدّت فرسایش درحالت کی را برای واحد مربوطه بیان می‌کند . برای دقیق‌تر شدن نتایج» سه پدیده 
فرسایشی سطحی (8). شیاری (18) و آبکندی ()) مد نظر قرار گرفتند. زیرنویس اوّل یکی از طبقات ۱ تا ۵ 
را شامل می‌شود و بیان‌گر شدّت فرسایش واقع بر این پدیده‌ها می‌باشد. زیرنویس دوّم حاصل سطح پوشش 
هر پدیده بوده و یکی از کدهای ۱ تا ۴ را به‌حود احتصاص‌داده و به‌ترتیب برای ۰ تا ۲۵ ۲۵ تا۵۰ ۵۰ تا ۷۵ 
و سطوح پوششی ۷۵ تا ۱۰۰ درصد استفاده شد. به‌طور کلی یک واحد کاری با کدی مشل 523121042 
نشان داده می‌شود. از سویی. برای اندازه‌گیری‌های صحرایی نیز با توجه به واحدهای کاری به‌دست آمده اقدام 
به نصب پلات‌های سطحی در هر واحدکاری با سه تکرار و به ابعاد یک متر و مساحت یک مترمرتع شد. در 
مورد شیارها و آبکندها از روش صادقی (۱۳۸۸) استفاده شد. در مورد شیارها نیز در هر واحدکاری باسه 
تکرار شیارهایی با طول بیش از ۰ متر انتخاب و مقطع عرضی آن‌ها در در سه مقطع ۲۵ ۰ و ۷۵ درصد از 
طول آن‌ها اندازه‌گیری‌شد. در مورد آبکندها از پیکه‌های چوبی استفاده‌شد تا دقّت مورد نظر تأمین گردد. در 
طول دوره یک‌ساله. دوبار اقدام به اندازه‌گیری رسوب شیارهاء آبکندها شد. همچنین در مورد پلات‌های 
سطحی نیز رواناب و رسوب حاصل از بارندگی‌ها به صورت تجمعی در چند نوبت مورد اندازهگیری قرار 
گرفت (شکل ۳ نقشه شدات فرسایش در واحدهای کاربری‌اراضی را نشان می‌دهد). 
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3 0.015 0 زباد #* 


سیر 


30515107۸ 30517070 30۶8190701 30۶2110701 30۶23103 0 ۳ 


ك 


4 :9010 کم 


...۰ وود 


1 


1 


1 


1 


1 


راهنما 
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مرتع. باغ و جنگل 


۴ بحث و نتایج 


نتایج آزمون آماری کروسکال والیس در بررسی توان ردیاب‌ها در جداسازی واحدهای کاری شامل 


اراضی مرتعیء اراضی جنگلی» اراضی زراعی و باغات در منطقه تنگ بستانک در جدول ۱ ارائه شده است. 
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شماره پیست و دوم 


سال شش استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت ... ۳۷ 


جدول ! نتایج آزمون کراسکال - والیس جهت تعیین خصوصیت منشأّیاب مناسب برای تفکیک منابع رسوب 


منشایاب آماره ۲ سطح معنی‌داری 
باریم ۶2۳۸/۸ ۰۱۳۸۹۱ 
کادمیوم ۱/9۸۳ ۰-۳۳ 
کبالت ۷/۵۹ ۰۳۹۸ 
کروم ۸۱/۹۷۴ ۰/۷۵ 
مس ۱۳/۷۳ ۰۳۴ 
یت 2۴۲ ۰/۴۸۸ 
منگنز 2/۸۵۵ ۴۰ 
نئودیمیوم ۹/۹۴۵ ۰۱۱۳۷ 
نیکل ۱/۳۳۶ ۰۶۸ 
۱۹۳ ۳ 
سیلیکون ۱۱5/۶۰۷ ۰/۰۳۸ 
استرانسبوم ۱۳/۴۴۵ ۰۴ 
تیتانیوم ۱۵/۰۹ ۲ 
وانادیوم 1/۴ ۳۸( 
ننک ۱/۷۴ ۰,9۶۷ 
نیتروژن تاه( ۰/۰۵ 
کزین ۱۶/۱۷۸ ۳ 
استرانسیوم ۸۷/ استرانسیوم ۸۶ ۳/۱ ۰۷۹۶ 
نئودیمیوم ۱۴۴/ نئودیمیوم ۱۴۳ 2۹۵ ۰۳۲ 


ازنظر این آزمون اعتلاف میانگین‌های هریک از عناصر درواحدهای کاربری‌اراضی به غیر از شش عنصر کربن» 
نیتروژن. مس. استرانسیوم. تیتانیوم و سیلیکون درسطح کمتر از ۵ درصد معنی‌دار نیست. در واقع تمام عناصر به 
غیر از شش عنصر یاد شده دارای قدرت تفکیک و جداسازی واحدهای مذکور نیستند. نتایج تحلیل تشخیص 
(0۳۸) و معنی‌داری آماری نشان داده ترکیب بهينةٌ خصوصیات منشأًیاب باید طوری انتخاب شود که علاوه بر 
داشتن کمترین هم‌حطی چندگانه» بیشترین توان تفکیک منابع را داشته باشد تا باعث مشکلاتی درتفسیر اهمیّت هر 
یک از متغیّرها و ناپایداری نتایج طبقه‌بندی تحلیل نشود. معیار کمترین هم خطی برای عناصر ترکیب‌ها داشتن 
عامل تورم واریانس حداقل ۱۰ است. ضرایب همبستگی بیشتر از ۸۹۵ با حذف یکی از دو ماتریس همبستگی 
اصلاح می‌شود. جدول (۲) مراحل مختلف اضافه شدن عناصر ردیاب در توابع تشخیص برای واحدهای کاربری 
اراضی به روش گام به گام و تأثیر آن‌ها بر توان جداسازی تحلیل تشخیص را نشان می‌دهد. با اضافه شدن هر 
عنصر مقدار مریّم فاصلهٌ ماهالائوبیس افزایش و میزان 1200002 ٩‏ ۷۷1 کاهش می‌یابد. چهار تابع بر اساس 
عناصر یاد شده درواحدهای کاربری اراضی استخراج شده‌اند و جدول (۳) معنی‌داری توابع را به وسیله آماره 
مریم کای نشان می‌دهند که توابع اوّل در واحدهای تفکیک منابع رسوب ۶۵/۳/ را به حود اختصاص داده‌اند. توابع 


۳۸ 


جغرافیا و مخاطرات محیطی 


شماره پیست و دوم 


او رت هت کی کاترن سر تیاو االاف‌ها را کرزها تفاس کل فومت طفنیی 


عامل ت 
ِ ور 
۳ خصوصیات منشایاب حداقل مربع فاصله 0 عدد تحمل 1 وزن ردیاب 
‌ واریانس 
وارد شده ماهالانویب ۵ (ععصهه۲0۱۵) منشایاب در من 
ِ ۳10 3 
۱ کرت 2۳5 ۱(کربن) ۱ ۱/۸ 
۸ کرین) ۱۶ ۱/۵ 
۲ ی ۳۹۱ ۳۵۱ # 
۸مس) ۱/۳۶ 
۶ کرین) ۱/۳ ۱/۳ 
۳ سپلیکون ال و ۸(مس) ۱/۳۳۶ 
۷ سیلیکون) ۱/۳۳ 
۴ (کرین) ۱۳۶ ۱/۵ 
۷ مس) ۱/۳ 
۴ تیتانیوم ۱/۳۳۳ ۰/۱۶۹ 
۸ سیلیکون) ۱/۲ 
۶ (نیتانیوم) ۱/۸۵ 
حدول۳ مشخصات توابع تشخیص حاصل برای تفکیک منابع رسوب (کاربری اراضی) 
مقدار درصد درصد تب همبستگ : آماره 
۳ تجمعی ۱ تا 
ویژه واریانس واریانس کانونی ۱2 مربع کای 
۱ ۱0۷ 2/۳ 2/۳ ۰/۷۸۲ ۰/۱۶۹ 2۳۸۳۶ ۰/۰ 
۲ ۰/۷ ۳۲ ۸۷/۳ ۳۶ ۳۰/۲۸ 1-۱ 
۳ ۰/۳۹۵ ۱/۳ ۹/۶ ۳۷۷ ۷۴۵ ۱۰۸/۵۸ ۰-۶ 
۴ ۰/۰۳۷ ۱/۴ ۱۰۰ ۰/۱۸۸ ۰/۹۶۵ ۱۳۴ ۰/۷۳۶ 
همان گونه که مشاهده شد تعدادی از عناصر در مرحله اول (۶عنصر) توان تفکیک منابع دانشان داد اه ان سمل 


نیتروژن و استرانسیوم بودند که این عناصر در مرحله دوم و در تابم تشخیص حذف شدند و در جمع ترکیب بهینه 
قرار نگرفتند. از مجموع ۱٩‏ عنصر و نسبت ایزوتوپی» تنها ۴ عنصر مس کربن. سیلیکون و تیتانیوم (عناصر نادر) قادر 
به جداسازی منابع رسوب بودند. بدین‌ترتیب نتایج تحلیل تشخیص نشان داد ۴ ردیاب آلی و ژئوشیمیایی قادر به 
تفکیک منابع رسوب در حوزه آبخیز تنگ بستانک بودند. ضرائب تصحیح مواد آلی در منابع رسوب و ضرائب وزنی- 
ویژه عناصر انتخابی و ضریب تصحیح سطحویژه ذرات در جدول ۴ و ۵ ارائه شده است. بابه حداقل رساندن 
معادلات ۴ ۸۵ ۶ ۷ و ۸ برای هر یک از نمونه‌های رسوب و متوسط‌گیری از مقادیر حاصل سهم هر یک از منابع 


لول عربه تست آمله اس همان‌گونه که مشاهده می‌شود در تمامی این جداول سهم کل بر حسب درصد در 


سال * 


استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... 


جدول۴ ضریب تصحیح کربن آلی (ز(6) و ضریب تصحیح سطح ویژه ذرات (ز) 


منابع تولید رسوب ضریب تصحیح کربن آلی 
باغات ۰۱-۹ 
کشاورزی 1/۶۸ 
جنگل ضعیف ۰۱۱۶ 
مرتع ۶ 


ضریب تصحیح سطح ویژه ذرات 
۰/۹۳۵ 


۰/۹۵۵ 
12 


۰۸۳۸ 


جدول ۵ ضرائب وزنی ویژه عناصر (:۷۷) کاربری اراضی 


ترکیب بهینه خصوصیات منشأیاب 


ضریب وزنی ویژه 


۰۳۵۷ 
۰۱۳۸ 


"۱۱۸۳۲ 
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جدول۶ سهم نسبی کاربری اراضی در تولید رسوب بر حسب درصد 


مدل کالینز (وصتلام)) 

۶2۷۷ 

۱/۲ 

۱۳/۳۸ 

۳۷/۹۳ 

مدل کالینز اصلاح‌شده (00[1105 ۷۲0011160) 

2۳۲ 
۹/۸۳۱ 

۱۷/۳۳ 

(2 

مدل لاندور (12800/607) 

۰/۲ 

۱۸/۵۹ 

۱۹/۵۲ 

*<:«./(۹۶ 

مدل موتا (۷۲012) 

۹/۰۸۰۲ 

۱/۹۷۵ 

۱/۸۰۵ 

۱۹/۶۳۶ 

مدل اسلاتری (51211677) 

۵/۶۸ 

۱۳/۳۲ 

۱۹/۶۴ 

2/۸۷ 


نوع کاربری اراضی 
زراعت ابی 
اراضی جنگلی 
اراضی باغی 
اراضی مرتعی 
نوع کاربری اراضی 
زراعت ابی 
اراضی جنگلی 
اراضی باغی 
اراضی مرتعی 
نوع کاربری اراضی 
زراعت آبی 
اراضی جنگلی 
اراضی باغی 
اراضی مرتعی 
نوع کاربری اراضی 
زراعت آبی 
اراضی جنگلی 
اراضی باغی 
اراضی مرتعی 
نوع کاربری اراضی 
زراعت آبی 
اراضی جنگلی 
اراضی باغی 
از ای کرت 


۳۰ جغرافیا و مخاطرات محیطی شماره پیست و دوم 
همان‌طور که از این جدول مشخص هست مدل کالینز اصلاح شده بیشترین مقدار 00 و کمترین مقدار 
۷۲ را به خود اختصاص داده‌است. بسیاری از محققان از جمله حدادجی و همکاران (۰)۲۰۱۵ کالینز و 
همکاران (۲۰۱۰ و ۲۰۱۲) وجود عامل زز۷٩‏ را عامل بهتر شدن این مدل نسبت به مدل کلاسیک و قبلی آن 
عنوان کرده‌اند و همان‌گونه که مشاهده می‌شود شاخحص 0600۳ و ۷15 نیز این مهم را نشان می‌دهد. 
مطابق با این روش بهترین مدل‌ها برای بیان سهم منابع رسوب مدل‌های کالینز اصلاح‌شده و بر اساس مقادیر 
اندازه گیری‌شده مدل موتا (۷۲012) هستند. جدول ۶ نتایج سهم منابع مختلف رسوب را برای واحدهای 
کاربری اراضی تیان می‌دهد. در تحقیقات داحلی بیشتر مدل کالینز کلاسیک و کمتر مدل کالینز اصلاح‌شده 
جهت محاسبه سهم منابع مختلف رسوب استفاده شده و به‌ندرت به کاربرد مدل‌های دیگر پرداخته شده است 
در حالی که نتایج نشان داد مدلی همچون موتا (۷۲۵92) نتایج نزدیک و قابل قبولی نسبت به مقادیر اندازه- 
گیری شلاه و نسبت به مدل پذیرفته شده کالینز اصلاح‌شده نشان می‌دهد. مدل کلاسیک کالینز که در بسیاری 
از تحقیقات داخلی ازران انتفاقه شده ات دفت کمتری دز مقایسه با نسخه اصلاح‌شده 1 از خود نشان داده 
کشاورزی و بیشترین سهم فرسایش و رسوب توسط مراتع می‌باشد. 


حدول ۷ نتایج حاصل از )0 و ۳ برای مقادیر بهینه‌سازی شده (1,0691) 


نام مدل برت) ۳۴9۲ 
موتا (۷۲0/2) > ۵۳۹ 
اسلاتری (1210617) > ۵۱۶۸ 
لاندور (1۵00/61) ۳۳/۷ ۴/۵۵ 
کالینز (011105)) ما 0( 
کالینزاصلاح‌شده (۷۲)00[11185) ۹9۹/۹ ۳/۳۰۷ 


نتایج داده‌های صحرایی و اندازه گیری میدانی برای انواع کاربری و در کدواحدهای مختلف و برای 
تکرارها و فرسایش‌های سطحی, شیاری و خندقی در جدول ٩‏ آورده شده است. با توجه به نتایج اندازه‌گیری 
میدانی و سطح هر کاربری میزان و سهم مشارکت هریک محاسبه و با توجه به درصد مساحت هر کاربری 
اراضی از سطح کل حوزه آبخیز اهمیت‌نسبی آن‌ها محاسبه شده‌است. جدول ۸ بیانگر این مهم 
می‌باشلد 


سال 


استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... 


جدول۸ سهم و اهمیت نسبی واحدهای کاریری اراضی حاصل از اندازه گیری صحرایی 


همچنین مقدار اهمیّت‌نسبی کاربری‌های زراعت. جنگل. باغ و مراتع بر اساس مدل کالینز اصلاح‌شده به- 


جدول٩‏ نتایج نسبی تولید رسوب حاصل از انواع فرسایش (تن در هکتار) 


نام کاربری اراضی درصد مساحت 
تراعت ۳/۴ 
باغات ۳/0 
جنگل 2/۸ 
مرتع ۱۶/۷۵ 
کاربری | میانگین فرسایش | میانگین فرسایش 
اراضی سطحی شیاری 
۰/۷۵ ۰/۹ 
زراعت 
۰/۰۹۷۵ ۰/۳۰۲۵ 
۲/۵ ۲/۳۹ 
جنگل ۲/۷۸ ۳/۵۲ 
۲/۷۸ ۳/۵ 
۱۹۶ ۳/۳۵ 
باغ ۳/۴ ۴/۸۵ 
۳/۳۲ 2/۷۵ 
۱/۰۸ ۱۳۸۳ 
مرتع ۱۸/۶۷ ۱۶2/۷ 
۱۹/۲ ۷/۶۶ 


درصد مشارکت در تولید رسوب و فرسایش 


۰/۸/۱ 
۳/۳۳ 
۱۳۱/۵۵ 


۶/۸۹ 


میانگین فرسایش 


۰,۳۷۶ 
۶۳ 
۰/۶۵ 

۰۳۵ 


۰/۹۵۵ 


کد 


واحدکاری 


۳ 1 
۳ 1 


۳ 2 
۳۵ 1 
]1( 1 


۳ 1 
۵ 07 
۱ 


(۱ 1 
۵ 2 
]1( 02 


آهمیّت‌نسبی 
۰/۰۳۱ 
۷۹۹ 
۰/۵۵۱ 


۳/۸۷ 


2/۵۱۹ 


1/۹۹۳ 


۳۸/۳ 


۳۱ 


ترتیب ۰/۲۸۲ ۰/۱۶۲ ۷/۷۲ و ۳/۴ بود. با توجه به نتایج اهمیّت‌نسبی حاصل از اندازه‌گیری میدانی سهم 


منابع تولید رسوب و نتایج حاصل از گزینش بهترین مدل ترکیبی نتیجه گیری کرد که اراضی مرتعی و باغی 


بیشترین اولویّت را نسبت به دیگر کاربری‌ها به حود احتصاص داده‌اند. مخصوصاً اراضی باغی که وسعت 


بسیار کمتری در مقایسه با اراضی مرتعی دارند و سهم بیشتری از فرسایش و رسوب را در مقایسه با کاربری 


مرتع روانه خروجی حوزه می‌کنند. 


۳ جغرافیا و مخاطرات محیطی شماره پیست و دوم 
۵ نتیجه گیری 

حضور عنصر کربن (/6) به‌عنوان ردیاب در تفکیک منابع رسوب (کاربری‌های‌اراضی) دیده شد. که با 
نتایج والینگ و همکاران (۲۰۰۸) کالینز و همکاران (۲۰۱۰), ژانگ و همکاران (۲۰۱۳), حددادچی و 
همکاران (۲۰۱۳)» ونگ و همکاران (۲۰۱۴) و نصرتی و همکاران (۱۳۹۴) مبنی بر حضور کربن به‌عنوان یک 
ردیاب ثابت در تفکیک کاربری‌های اراضی مطابقت می‌کند. حضور عنصر کربن آلی در تمایز بین کاربری- 
های مختلف می‌تواند نشان‌دهنده تأثیر حاصل خیزکننده‌ها باشد. نتایجم حکیم‌خانی و همکاران (۱۳۸۷) مبین 
نقش متمایزکننده عناصر نادرخاکی از جمله سیلیکون. تیتانیوم و مس در تفکیک منابع رسوب می‌باشد. عدم 
حضور عنصر نیتروژن در مرحله دوّم انتخاب بهترین ترکیب بهینه ردياب‌ها با استفاده از تابع تحلیل تشخیص 
گام به گام می‌تواند ناشی از هم‌خطی بین ردیاب‌ها باشد که نتایج تحقیقات والینگ (۲۰۰۵) و نوحه‌گر و 
همکاران (۱۳۹۵) به این مهم اذعان دارند. بیشترین اهمیّت‌نسبی مربوط به کاربری مراتع و سپس باغات می- 
باشد. کاربری زراعت به دلیل قرار گرفتن در دشت بسیار مسطح (شیب بسیار کم» کشت متناوب (برنج و 
گندم) و بدون آیش در منطقه دارای کمترین اهمیّت‌نسبی از لحاظ تولید رسوب در محدوده مطالعاتی می- 
باشند. شاحص اهمیّت‌نسبی سهم‌های فرسایشی اولویّت‌بندی مدیریت کاربری اراضی را نشان می‌دهده 
چنان‌ که مراتع بالاترین اولویّت را در بحث کنترل تولید رسوب و لزوم مدیریت صحیح دارند. همچنین نتایح 
حاصل از تحقیق حاضر حکایت از بیشترین دقت در تبیین سهم‌های فرسایش توسط مدل کالینز اصلاح‌شده 
داشت که با نتایج حدادچی و همکاران (۲۰۱۵) مبنی بر انتخاب این مدل به‌عنوان یکی از قوی‌ترین مدل‌های 
ترکیبی در تعیین سهم منابع مختلف رسوب مطابقت دارد. تفاوت بین مقادیر اندازه‌گیری شده و محاسباتی 
مدل‌های ترکیبی (روش انگشت‌نگاری) ناشی از تله‌اندازی رسوب و نسبت تحویل آن می‌باشد که سبب 
تفاوت در مقادیر این دو روش می‌شود. هرچند هر دو روش از حیث روند همدیگر را تأْیید می‌کنند. تفاوت - 
های مقداری حاصل از دو روش را می‌توان به مباحثی چون نسبت تحویل رسوب کلّی و به تفکیک منابع 
رسوب. ذخیره جابجایی و انتقال رسوب نسبت داد. همچنین تفاوت بین سهم‌های مشاهداتی و محاسباتی با 
مدل‌های کمّی را می‌توان به ضعف روش‌بهینه‌سازی لوکال نسبت داد که فقط با در نظر گرفتن جواب‌های 
محلّی اقدام به بهینه‌سازی تابع هدف می‌کند و مانند دیگر روش‌های بهینه‌سازی از جمله الگوریتمژنتیک 
(0۸) کلیّه جواب‌های عمومی و انتخاب بهترین جواب را مدنظر قرار نمی‌دهد؛ بنابراین پيشنهاد می‌گردد در 
مطالعات منشأٌیابی و بهینه‌سازی مدل‌های ترکیبی از دیگر روش‌های بهینه‌سازی که دقت بیشتری دارند 
استفاده شود. اگرچه در این تحقیق کارآیی مدل‌های ترکیبی با روند اندازه‌گیری صحرایی مطابفت داشت 
(مدل موتا با ضریب همبستگی ۰/۹۲۴ مدل اسلاتری با ضریب همبستگی ۰/۸۹۲ مدل کالینز با ضریب 
همبستگی ۰/۷۶ مدل لاندور با ضریب همبستگی ۰/۸۳۴ مدل کالینز اصلاح‌شده با ضریب همبستگی 


بالیشتم استفاده از مدل‌های ترکیبی و ردیاب‌ها جهت... ۳۷ 
۳ با در نظر گرفتن دو شاخص 00 و ]۷[ نتایج مدل کالینزاصلاح‌شده با وجود ضریب 
همبستگی کمتر مقبول می‌باشد و نتایج سهم‌بندی مدل موتا (0۷10/2) و با اسلاتری (۹1211610) با وجود 
شاخص منفی 00 غیرقابل قبول هست (هر چه مقادیر شاخص 0۳ به یک نزدیک‌تر باشد نتایج مدل 
قوی‌تر و فابل اعتماد می‌باشد). اگرچه این مدل توانسته ضریب همبستگی بالاتری با مقادیر اندازه گیری شده 
نشان دهد. یکی از موارد ضعف مدل‌های ترکیبی در نشان دادن سهم منابع تولید رسوب عدم در نظر گرفتن 
مواردی همچون سطح منابع می‌باشد. برای رفع این مهم از شاخص اهمیّت‌نسبی جهت مقایسه منابع مود 
رسوب و مقایسه یکسان آن‌ها استفاده می‌شود. در بیشتر تحقیقات داخلی مقایسه‌ای بین مدل‌های کمّی 
مختلف صورت نگرفتهاست و پا در تحقیقات خارجی بجز چندمورد مقایسه‌ای بین نتایج این مدل‌ها انجام 
نشده‌است. در این تحقیقات جهت گزینش بهترین مدل از شاخحص 26(۳) استفاده‌شده و کمتر به شاحص‌های 
دیگر از جمله ۷15[ توجه شده‌است؛ بنابراین استفاده از این شاخص‌ها جهت تأیید نتایج و نیز مقایسه 
این نتایج با مقادیر اندازه‌گیری شده ضروریست. 


کتاینامه 
کارشناسی‌ارشد دانشگاه زابل دانشکده آب و خاک گروه مرتع و آبخیزداری. 
صادقی» س.ح.ر. (۱۳۸۸). مطالعه و اندازهگیری فرسای شسآبی, انتشارات دانشگاه تربیت مدرس. 
حکیم‌خانی؛ ش. (۱۳۸۹). ارزیابی اهمیت نسبی انواع فرسایش در تولید رسوبت (بررسی موردی حوزه قره‌اقاج ماکو). 
حیدری» ک.. نجفی‌نژاده ع.» خرمالی؛ ف. و بابانژاد. م (۳۹۲). تعیین سم واحدهای کاری در تولید رسوب معا 
استفاده از روش منشاأیابی رسوب (مطالعه موردی آبخیز تول‌بنه, استان گلستان), فصلنامه پژوهش‌های فرسایش 
محیطی. ۳ (۱۱» صفحات ۲۷ تا ۳۸ 
عباسی. 0 فیض نیا س-۰ احمدی. ح و کاظمی. ی (۱۳۸۵). منشایابی رسویات بادی نیانک با استفاده از ردیابی 
غلامی» ح.» طاهری‌مقدم» قیری. م. و مهدوی نجف‌آبادی» ر. (۱۳۹۴). تعیین سهم کاربری‌های اراضی در تولید 
رسوب تبه‌های ماسه‌ای با استفاده از روش انگشت‌نگاری (مطالعه موردی: ارگ نگار بردسیر کرمان) مجله 
سیمای فرسایش» لکشت تکاوضن و اندازه‌گیری میدانی» محله مهن سی و مدیریت آبخیز ۵ (۳) صفحات ۱۶۵ تا 


۷/۸ 


۳۴ جغرافیا و مخاطرات محیطی شماره پیست و دوم 


معلق و کف بر پایه منشأیابی رسوب در حوزه طالقانی» مجله منابع طبیعی ایران ۶۸ (۴ صفحات ۷۵۱ تا ۷۶۵. 
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